




















OF MICHiCan 
MAY 07 igo) 





The el 
Vieteorological 


Magazine 


March 1991 
Vol. 120 No. 1424 





551.577.21(495) 





Spatial distribution of rainfall in the Greater 


Athens Area 


G.T. Amanatidis, C. Housiadas and J.G. Bartzis 
National Center for Scientific Research Demokritos, Institute of Nuclear 
Technology and Radiation Protection, 15310 Ag. Paraskevi Attikis, Greece 





Summary 


The annual and seasonal spatial distribution of rainfall in the Greater Athens Area in Greece is studied using the 
rainfall data from 24 rain-gauges for 1985-89. Both these rainfall distributions and the temporal variations of 
rainfall during the 5-year period present strong discrepancies among the different parts of the examined area. The 
influence of the complex topography and the atmospheric circulation on the rainfall distribution pattern for 
different seasons is examined. Finally, the relation of rainfall amounts to the station height above mean sea level is 


also investigated. 


1. Introduction 

The knowledge of the spatial distribution of rainfall is 
of interest not only from a meteorological viewpoint, 
but also for its importance in different fields such as 
agriculture, hydrology, water resources, atmospheric 
pollution or even in flood control. The estimation of 
the spatial distribution of rainfall is a complex task, 
especially in cases where detailed information concern- 
ing the impact of local topography on the dominant 
atmospheric circulation is not available. The problem 
becomes even more complicated by the varying pattern 
and intensity of near-surface atmospheric circulations. 

In the case of the Greater Athens Area (GAA) (see 
Fig. 1) the knowledge of the spatial distribution of 
rainfall is rather limited though rainfall measurements 
have been carried out in the city of Athens for well over 
130 years. One of the longest series of rainfall 
measurements in south-eastern Europe is that collected 
by the Meteorological Institute of the National Observ- 
atory of Athens (MINOA). This has been extensively 
analysed in the last 15 years (Katsoulis et al. 1976, 


Zerefos et al. 1977, Repapis 1986, Katsoulis and 
Kambezidis 1989), but since the knowledge is based on 
rainfall data obtained from only one rain-gauge, 
nothing can be said about the spatial distribution of 
rainfall in the GAA. Extrapolation of the results 
obtained at MINOA could create serious errors caused 
by the complex topography effects, especially during the 
cold period of the year when frontal rainfall is more 
common. During the warm period of the year, rainfall is 
predominantly convective and frequently has high 
intensity over small areas. Such local thunderstorm 
activity makes data extrapolation from one station even 
more unreliable. 

In this study, the rainfall measured by 24 rain-gauges 
in the GAA for the 5-year period 1985-89 are analysed, 
in order to study the spatial distribution of rainfall in 
this area. The network used consists of all the reliable 
rain-gauges in operation for the examined period and 
covers the area in a mesoscale range. Although the 
results do not have an adequate climatological signific- 
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Figure 1. Map of the Attika peninsula. The Greater Athens Area is 
indicated by the shaded domain. 














ance, because of the relatively short examined period, 
they confirm the expected strong irregularity of rainfall 
distribution in the GAA. 

The results of this study may be useful in the 
assessment of numerical models that aim at predicting 
the rainfall pattern in the GAA. One such model is 
currently under development at the National Centre for 
Scientific Research (NCSR) ‘Demokritos’ for simulating 
rainfall distribution in complex topography (Housiadas 
et al. 1991). This model constitutes a further improve- 
ment of the three-dimensional atmospheric transport 
code, ADREA-I, developed in the context of environ- 
mental impact assessments for wind field and dispersion 
calculations (Bartzis et al. 1990, 1991). 


2. Topography 

The area studied in this paper is located at the south- 
eastern part of the Attika peninsula. The presence of 
mountains and sea in the GAA divide it into different 
regions; namely, the Athens Basin, the Thriassio Field, 
the Mesogia Plain, the Marathon Area and the islands 
of Salamina and Egina (Fig. 2). 

The Athens Basin, where the city of Athens is located, 
is surrounded by mountains on three sides, while to the 
south there is sea (Saronikos Gulf). As shown in Fig. 2, 
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Figure 2. Map of the Greater Athens Area. Land above 400 m is indicated by shading, while the rain-gauge stations (see Table I) are indicated 


by numbered stars. 
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Mt Aegaleo (468 m) is situated to the west, Mt Parnitha 
(1413 m) and Penteli (1107 m) to the north and Mt 
Hymettos (1026 m) to the east. The Athens Basin is 
30 km long and 12 km wide with a south-west to north- 
east major axis, and is interspersed by a series of small 
hills. 

The Thriassio Field, located at the west side of the 
Athens Basin, is also surrounded by mountains. Mt 
Aegaleo is situated to the east, Mt Parnitha to the north, 
Mt Kitheron (1409 m) and Mt Pateras (1131 m) to the 
west, while to the south is the Elefsis Gulf coastline. 
Inside this area most of the Greek industry is 
concentrated, mainly along the coast, while inland, 
agricultural activities are dominant. To the east of the 
Thriassio Field and south of Mt Pateras extends the 
Megara Area, part of which is also included in the 
GAA. 

The Mesogia Plain is located at the east of the Athens 
Basin and is separated from it by Mt Hymettos. 
Mt Penteli is situated to the north, Mt Panio (636 m) to 
the south, while to the east there is the Aegean Sea. 
Inside the plain there are some sparse industries, but 
agricultural activities predominate. 

The Marathon Area is located to the north of Mt 
Penteli and east of Mt Parnitha and is washed to the 
north and east by the Evoikos Gulf. In the centre of the 
area there are the Marathon dam and the artificial 
homonymous lake. The lake is 250 m above mean sea 
level and is surrounded by forested hills. 

Finally, the islands of Salamina and Egina are 
situated at the south-eastern part of the GAA in the 
Elefsis and Saronikos Gulfs respectively. 


3. Climate 

The climate of the GAA is Mediterranean with hot, 
dry summers and wet, mild winters. The summer dry 
period lasts for about 6 months. It begins at the end of 
April and continues until the beginning of October. The 
average daily winter temperature is 9.9°C and the 
summer one is 25.8 °C. The total solar radiation is rather 
strong with average daily values of about 22 MJ m” in 
the summer and 8 MJ m” in the winter. 

Most of the annual rainfall occurs in the months from 
October to March. This is due to the synoptic 
circulation which is related to extensive disturbances 
within the zonal westerlies (Repapis 1986). More 
precisely, during the summer the Atlantic anticyclone 
with its more northerly position dominates the Medit- 
erranean region and therefore deflects the tracks of the 
storm systems. Consequently, summer is an inactive 
season and only infrequent local thunderstorms occur, 
especially in the afternoon. During the cold period of the 
year, this high moves to more south-westerly positions 
allowing depressions to enter the Mediterranean or to 
develop within the region. These rain-bearing depressions 
move from west to east (Conrad 1943, Meteorological 
Office 1962) and are responsible for great rainfall 
amounts over Greece. Contributions to the rainfall over 


Greece are also provided by the extension of the 
Siberian anticyclone over the Balkans. This creates a 
cold northerly airstream which, combined with possible 
warm and wet air masses over the central Mediterranean, 
produces considerable amounts of rainfall. 

The upper-air flow over Greece is characterized by 
prevailing wind directions from the west/south-westerly 
sector during winter and spring, and of the northerly 
sector during summer and autumn (Pissimanis et al. 
1989). Over the Athens Basin the winds tend to blow 
along the axis of the basin (Andreakos et al. 1984). The 
most frequent wind is that blowing from the northerly 
sector. The second most frequent wind blows from the 
south sector, while those from the east and the west 
sectors have a very low frequency. Except for the 
synoptically driven north or south winds, there are two 
permanent weather systems of regional character that 
substantially contribute to this high frequency of 
northerly and southerly winds. Firstly, the combination 
of high barometric pressure over the Balkans with low 
pressure over the eastern Mediterranean causes the 
north winds, known as the ‘Etesians’, in the summer 
months (Meteorological Office 1965). Secondly, the 
local sea- and land-breeze mechanism is responsible for 
winds from the southern sector during the warm hours 
of the day in late spring and in summer (Prezerakos 
1986). Another sea-breeze cell is established in the 
Mesogia Plain and gives rise to an airflow through the 
pass between the Hymettos and Penteli mountains. At 
that point, because of the vertical motion created by the 
opposing sea-breeze currents (south-west and north- 
east), some convective clouds appear that are occasion- 
ally significant (Clement 1989). 


4. Rainfall data in the Greater Athens Area 

The 24 rain-gauges selected for the present investig- 
ation belong to different Greek institutions: the 
previously mentioned MINOA station, ten stations of 
the Water and Sewage Corporation (WSC), six stations 
of the National Meteorological Service (NMS), five 
stations of the Ministry of Environment, Planning and 
Public Works (MEPPW), one station of the Vineyard 
Institute of the Ministry of Agriculture (VIMA) and one 
private station. The list of the stations used, together 
with their geographical coordinates and their elevations 
above mean sea level, is given in Table I, and their 
network distribution is shown in Fig. 2. 

Since the objective of the present study was to 
determine the spatial distribution of rainfall the GAA, a 
major effort was made to take into account data from 
the maximum number of stations with continuous 
records for the period 1985-89. On the other hand, 
special care was given to avoid effects arising from 
improper exposure, change of exposure, instrumental 
failures or change of station site during the examined 
period. The representativeness of the 5-year period was 
also examined by comparing it with the climatological 
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Table I. 


Rainfall stations in the Greater Athens Area used in this study. See Fig. 2 


for locations, and text for explanation of acronyms. 








No. Station Institution Longitude Latitude Height 
(m) 
1 MINOA MINOA 23° 43’E 37° 58’N 107 
2 Filadelfia NMS 23° 44’E 38° 02’N 136 
3 Peireas NMS 23° XE 37° 56’N 3 
4 Helliniko NMS 23° 44’E 37° 54’°N 10 
5 Eleysina NMS 23° 33’E 38° 03’N 30 
6 Megara NMS 23° 2VE 38° 00’N 36 
7 Egina NMS 23° 26’E 37° 45’'N 3 
8 Vyronas MEPPW 23° 45’E 37° STN 120 
9 Nikaia MEPPW 23° 39'E 37° 58’N 20 
10 Halandri MEPPW 23° 48’E 38° 02’N 180 
11 Peristeri MEPPW 23° 42’E 38° O1’'N 60 
12 Markopoulo MEPPW 23° 56’E 37° 53’N 85 
13 Anavyssos Private 23° S7’E 37° 42’N 80 
14 Lykovrysi VIMA 23° 47’E 38° 04’N 220 
15 Tatoi WSC 23° 47’E 38° 06’N 237 
16 Horio WSC 23° 58’E 38° 09’N 60 
17 Fragma WSC 23° 54’E 38° 10’N 240 
18 Stamata WSC 23" Se 38° 08’N 370 
19 Kapandriti WSC 23° 53’E 38° 13’N 340 
20 Varnavas WSC 23° SSE 38° 14’N 480 
21 Kiourka WSC 23° 50’E 38° 12’N 400 
22 Katsimidi WSC 23° 49’E 38° 11’N 570 
23 Hani WSC as SNE 38° 08’N 260 
24 Bogiati WSC 23° SIE 38° 08’N 350 





data obtained from the longest rainfall time-series of 
MINOA. 

In order to study the rainfall amounts in the different 
regions in the GAA (as defined in section 2), the rainfall 
data from the 24 stations were grouped into six 
categories for further analysis; the criterion for this 
grouping was mainly geographical. Because of the 
length of the Athens Basin, it was subdivided into the 
South and North Athens Basins. The height of each 
station above mean sea level was also taken into 
consideration in this division, as a relationship between 
station height and rainfall amount exists (Smith 1979). 

The categories into which the GAA was divided were 
the following: 

South Athens Basin (SAB) which includes six 

stations (MINOA, Peireas, Helliniko, Nikaia, Vyronas 

and Peristeri). 

North Athens Basin (NAB) which includes four 

stations (Filadelfia, Halandri, Lykovrysi and Tatoi). 

Thriassio Field which includes two stations (Eleysina 

and Megara). 

Egina Island with its homonymous station. 

Mesogia Plain with two stations (Markopoulo and 

Anavyssos). 

Marathon Area which includes nine stations (Horio, 

Fragma, Stamata, Kapandriti, Varnavas, Kiourka, 

Katsimidi, Hani and Bogiati). 


5. Results and discussion 


5.1 Temporal variations of rainfall 

The monthly values of rainfall in the six regions of the 
GAA for the period 1985-89 are shown in Fig. 3 and the 
corresponding monthly mean values for the same period 
are shown in Fig. 4. It is obvious, that the rainfall in the 
GAA presents a strong seasonal variation with wet 
winters and dry summers, as in the typical Mediterranean 
climate, which is a result mainly caused by the synoptic 
circulation explained before. The monthly values of 
rainfall vary greatly from month to month. Calculations 
from the long-term rainfall series of MINOA have 
shown that the correlation between pairs of months is 
very low, which has been attributed to the relative in- 
dependence among the synoptic conditions characterizing 
the rainfall in each season (Katsoulis and Kambezidis 
1989). It is important to note that 1989 was relatively 
dry, especially for the southern part of the GAA. In that 
year rainfall (150.6 mm) was the lowest measured at the 
MINOA station for this century. This effect is related to 
the general drought problem experienced recently 
(continuing during 1990) all over Greece, which has in 
particular created severe problems for the water 
resources of the city of Athens. 

The correlation among the records from the different 
parts of the GAA is relatively high, and the corresponding 
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Figure 3. The monthly values of rainfall in the six regions of the Greater Athens Area for the period 1985-89. 


coefficients are presented in Table II. The higher 
correlation coefficients are observed between the pairs 
of the SAB, NAB and Marathon Areas. This fact is 
explained by the coincidence of the geographically 
consecutive areas with the south-west to north-east 
alignment of the atmospheric circulation. The smaller 
correlation coefficient is observed when the area of 
Egina Island is considered. This may be attributed to the 
different behaviour of rainfall due to the airflow’s 


divergence—convergence mechanism commonly observed 
in island regions. 

The monthly mean variation of rainfall for the period 
1985-89 in the six parts of the GAA, presented in Fig. 4, 
is characterized by a persistent and unexpected maximum 
in March for all groups of stations. At the end of the year 
the maximum rainfall appears in November. Consequently, 
the 5-year period presents a secondary minimum during 
winter. Katsoulis and Kambezidis (1989) have also 
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Table Il. 


Area. See Fig. 2 for locations of areas. 


Correlation coefficients among the rainfall records in the Greater Athens 
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Figure 4. The seasonal variation of rainfall in the six regions of the Greater Athens Area. 


observed during the most recent period (1951-85) at the 
MINOA station an increasing tendency of rainfall 
amount in the transitional seasons (which are the 
unstable part of the year), although at the same time a 
noteworthy significantly negative trend was found to 
characterize the annual rainfall amounts. They qualified 
this phenomenon as indicative of an urban effect 
(Changnon 1968, Goldreich and Manes 1979). The 
slightly increasing tendencies of rainfall amount are 
attributed to the urban heat-island effect and the surface 
roughness in the city, which both increase the turbulence 
and the atmospheric instability. A similar conclusion 
was also drawn by Goldreich and Manes (1979) for the 


Tel-Aviv urban area. In this respect, it is interesting to 
note that in the GAA the less pronounced decreases 
appear at the Marathon and Egina Areas which are the 
less urbanized. 

In order to get more insight into the shift of the 
rainfall maximum towards March, the climatological 
data at the MINOA station for the periods 1891-1985, 
1951-80 and 1980-89 along with the period of interest 
1985-89 are plotted together in Fig. 5. The comparison 
shows that the March particularity appears also in the 
monthly mean variation of the last decade (1980-89). 
This result is likely to be of climatological interest. 
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Figure 5. 


5.2 Rainfall variation with station height 

As is well known, rainfall tends to increase with height 
and decreases with depth down in an isolated valley 
(Smith 1979). On very high mountains the rainfall 
increases up to a certain height (up to 2 km) and then 
decreases, because above this height the ascending air 
masses lose their water vapour content. The relationship 
between the average annual rainfall and the height of the 
station above mean sea level is in most cases linear. The 
scattering in the linear-regression representation is 
caused by the failure of the relationship to describe the 
upslope-rain to rain-shadow contrast, which is important 
for broad mountains. If the region under consideration 
has a consistent prevailing wind, this latter effect is well 
represented by the spatial distribution of annual 
rainfall. 

In Table III, the average annual rainfall and the mean 
height of the stations for the six different areas in the 
GAA are presented. It is ascertained that rainfall 
increases with height. The maximum rainfall occurs in 
the hilly area around the Marathon lake. On the other 
hand, the minimum rainfall occurs in the south-west 
part of the GAA at the Thriassio Field and the Egina 
Island where the stations are very close to sea level. 

In Fig. 6, the average annual rainfalls for the period 
1985-89 at the 24 stations in the GAA have been plotted 
against the station heights. Although the station height 
used is not obtained by any averaging process con- 
sidering the surrounding area (Goh and Lockwood 


The seasonal variations of rainfall at the MINOA station (see Table I) for the periods 1891-1985, 1951-80, 1980-89 and 1985-89. 





Table Ill. Rainfall amounts at different regions of the 
Greater Athens Area. The areas are the same as in Table II. 








Area Number of Average Average annual 
stations height rainfall 
(m) (mm) 
SAB 6 53 320.9 
NAB 4 193 465.8 
Mesogia 2 83 366.4 
Thriassio 2 33 275.5 
Egina 1 3 288.4 
Marathon 9 341 557.2 





1974, Hill et al. 1981) the correlation is relatively high 
(r = 0.826), which indicates that variations of annual 
rainfall amounts account for 68% of the variance of 
the station height. The 32% uncorrelated variance of 
rainfall is attributed to the upslope-rain to rain-shadow 
contrast, as mentioned before. The linear-regression 
equation, obtained in Fig.6, R= 306.8+0.685H is 
valid for the GAA, where the rainfall R is measured in 
millimetres and the station height H in metres. The slope 
of the linear best fit suggests that annual rainfall 
amounts increase by 68.5 mm for each 100 m increase of 
the station height. The intercept of the linear best fit 
indicates that there is an annual rainfall of about 
300 mm in the GAA for any stations with zero height. 
The correlation between the monthly rainfall amounts 
and the station height presents a seasonal variation. The 
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Figure 6. The average annual rainfall for the period 1985-89 versus 
station heights above mean sea level. 


correlation coefficients are high for January (0.696) and 
October (0.789), but are low for April (0.271) and July 
(0.155). These results confirm that the frontal rainfall 
during the cold period of the year is strongly affected by 
the presence of mountains, while the thunderstorm 
activity during the warm period of the year is less 
affected by the local topography. The October maximum 
correlation coefficient calculated for the GAA is in 
accordance with Bergeron (1968, 1973). During the fall 
months, he found in Uppsala, Sweden, a remarkably 
strong dependence of rainfall on height — primarily 
caused by rain from stratified clouds. 


5.3 Spatial variations of rainfall 

The great spatial variability characterizing the rainfall 
in the GAA is not only attributed to the synoptic 
disturbances. Geographical and topographical factors, 
such as mountains, hills, plains, sea and islands in the 
GAA, produce microclimatic effects on the mesoscale 
circulations and significantly disturb the distribution of 
rainfall in the GAA. As a result, there will be as many 
different spatial rainfall distributions as there are days 
considered. Some degree of generalization is therefore 
essential, and consequently the spatial distribution of 
rainfall in the GAA is analysed on a monthly basis. 

The average annual distribution of rainfall in the 
GAA for the 5-year period 1985-89 is shown in Fig. 7. 
The maximum rainfall occurs in the Marathon Area and 
more precisely on the south-east side of Mt Parnitha, 
while the minimum rainfall occurs in the south-west part 
of the GAA over the Salamina and Egina islands. A 
gradient is observed between inland and coastal areas 
along the south-west to north-east direction. The annual 
rainfall amounts double in the Marathon Area in 
comparison to the coastal areas. Repapis (personal 
communication) has also observed a similar gradient 
using earlier data from stations at SAB, NAB, Mesogia 
and Thriassio Areas. The rainfall gradient in the GAA 
is mainly attributed to the rain-bearing depressions 














Figure 7. The average annual spatial distribution of rainfall in the 
Greater Athens Area for the period 1985-89. 


moving from south-west to north-east over the Attika 
peninsula. The air masses are influenced by the presence 
of the surrounding sea which is warmer than the 
adjoining land during the cold period of the year. The 
passage from the sea to the land, together with the 
forced ascent due to topographic effects, are thought to 
be responsible for the gradient observed. Furthermore, 
superimposed on the above mechanisms, an additional 
effect seems to occur over the Marathon Area which 
explains the enhancement of the rainfall amounts 
measured in this area. Indeed, the hills and the small- 
scale mountains around the Marathon lake enhance 
rainfall by the so-called ‘seeder-feeder’ mechanism 
advanced by Bergeron (1965), in which raindrops from 
pre-existing (seeder) clouds aloft wash out small cloud 
droplets within low-level (feeder) clouds over the hills. 
According to field observations and theoretical studies 
(Bergeron 1968, 1973, Hill et a/. 1981, Carruthers and 
Choularton 1983, Grabowski 1989), hills receive 40— 
100% more rain than the surrounding areas, while for 
convective-type rainfall, they appear to have negligible 
influence on the spatial rain distribution. 

In order to analyse the seasonal aerial distribution 
of rainfall in the GAA, four months (January, April, 
July and October) considered as representative for the 
respective seasons are used. The spatial distribution 
for each month is shown in Fig. 8. The mean rainfall 
distribution for January and October present the same 
patterns as the annual one. But, in the January case, 
the gradient increases and the difference between the 
coastal area and the Marathon Area becomes more than 
three times bigger. During April, the rainfall pattern 
changes and the observed gradient becomes much less 
pronounced. Maximum amounts of rain fall on the 
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Figure 8. The spatial distribution of mean rainfall for the period 1985-89 in the Greater Athens Area for (a) January, (b) April, (c) July, and 


(d) October. 


foot of Mt Parnitha and minimum at the south side of 
Mt Hymettos. 

The July spatial distribution of rainfall is character- 
ized by very low values (dry period of year). The 
influence of topography on rainfall distribution decreases, 
since for this season rainfall is only due to the local 
thunderstorm activity. High values of rainfall are 
observed in the Marathon Area and at Helliniko. The 
latter is in accordance with the work by Prezerakos 
(1989), who observed a greater number of thunderstorm 
days at Helliniko. He claimed that the thunderstorms 
start over the continental parts of Athens Basin and then 
develop towards the coasts of Saronikos Gulf. 


6. Conclusions 
The spatial distribution of rainfall in the GAA in 
Greece was examined in relation to the complex 


topography and the atmospheric circulation of this 
area. For this purpose, rainfall data from 24 stations for 
the 5-year period 1985-89 was analysed. This analysis 
revealed the following features: 
(a) Although the rainfall in the GAA presents in 
general a seasonal variation with the maximum 
during the cold period of the year and a minimum 
during the warm period of the year, a secondary 
minimum is observed during winter in the examined 
period (appeared also in the last decade) intensified 
by the higher rainfall amounts measured in March. 
(b) The annual rainfall amounts increase with height 
in the GAA by 68.5 mm for each 100 m increase of the 
station height above the mean sea level. A background 
annual rainfall of about 300 mm is considered as 
independent from the relief in the GAA. The 
correlation coefficient between the monthly mean 
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rainfall and the station height presents a seasonal 
variation with high values during the cold period of 
the year and low values during the warm period. 

(c) The annual spatial distribution of rainfall presents 
a gradient between inland and coastal areas. As a 
result, rainfall at Marathon Area doubles in compari- 
son with the coastal areas. The same distribution 
pattern is observed for January and October, while 
for April and July the gradient is much less 
pronounced due to the thunderstorm activity, which 
is less influenced by the local topography. 
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A comparison of UK road ice prediction models 
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Summary 


The performance of three models are compared and analyzed using a variety of statistical techniques. 


1. Introduction 

The climate of the British Isles in winter is usually 
characterized by the fluctuation of road surface temper- 
atures around 0 °C together with a high humidity, which 
can induce the formation of ice or frost on road surfaces. 
This is obviously a serious potential hazard to alli road 
users. There is a need for an accurate prediction of road 
surface conditions, to allow time for winter maintenance 
engineers to spread de-icing chemicals to prevent the 
formation of ice or frost or the accumulation of snow. 
Such phenomena also occur in many other countries 
with acold and wet winter climate, such as northern and 
western Europe, Japan and North America. Three 
numerical road ice prediction models have been 
developed in the United Kingdom (Thornes 1984, 
Parmenter and Thornes 1986; Rayer 1987, Thompson 
1988; Shao 1990) specifically for this purpose. Some 
others have also been developed outside the United 
Kingdom, such as the Swedish (Kempe 1990), French 
(Isaka et al. 1990), American SSI model (Stephenson 
1988) and Finnish (Nysten 1980) models. Partly for 
commercial reasons, these other models are not 
available for comparison. 

There have been few comparisons between the road 
ice prediction models using observed road surface 
temperature data (Thornes 1989), the COST-309 (1991) 
study merely comments on the availability of such 
models. This paper presents the results of a detailed 


comparison between the three UK road ice prediction 
models and also suggests a standard set of statistical 
parameters for such comparisons using: 

(a) bias, standard deviation (SD) and r.m.s. error for 

overall, maximum and minimum temperatures, 

(b) freezing time and duration, and 

(c) forecast versus frost/no frost analysis. 


Table I shows, for a matrix of values of forecast minus 
actual road surface temperature, the different measures 
of bias, r.m.s. error and SD that can be calculated for 
hourly and daily values. 

For comparison, Chapman’s Hill road weather 
outstation on the MS motorway, West Midlands was 
chosen as the test site, due to the availability of roadside 
data and weather data from the University of Birmingham 
some 7 miles to the north-east. The three available 
models, i.e. the Thornes model (Thornes 1984), the Met. 
Office model (Rayer 1987, Thompson 1988) and the 
Icebreak model (Shao 1990) — developed on a 1987/88 
winter database — are compared using a standard set of 
input data for the site taken from the winter of 1988/89. 
A comparison of the similarities and differences of the 
three models is given in Table II. The models also use 
different methods for the computation of the radiative 
and turbulent fluxes. 





Table I. 


The derivation of hourly (H), daily (D) and overall (O), bias (B), r.m.s. error 


(R) and standard deviation (SD) for a matrix of values of forecast minus actual road 
surface temperatures. 





Hour of the day from noon (i) 





1 2 I 

Day (n) . 
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Table Il. 


Brief description of similarities and differences between the Icebreak, Thornes and Met. Office models 





Thornes Met. Office 





Icebreak 
1. Basic equation Heat conduction 
2. Road temperature Yes 
profile 
3. Energy balance Yes 
4. Longitude Yes 
5. Sky view factor Yes 
6. Influence of traffic Yes 
7. Inputs* T., Ta, W, Ca, C, R 
8. Outputs** T; and wetness 
9. Method of Control volume 
discreterization 
10. Derivation of Solving basic equation 
forecast by fully implicit 


scheme 


Energy balance 


Heat conduction 


No Yes 

Yes Yes 

Yes No 

No No 

No No 

Ty ta. Ge GR Te a ty ty Cy RP 
T; and wetness T; and wetness 


Searching for root of 
basic equation by 


Finite difference 


Solving basic equation 
by explicit scheme 


iteration 





* 


T, = air temperature; Ra = dew-point; W= wind speed; C, = total cloud amount; C, = dominant cloud type; C; = low cloud 


amount; R= rain period; P;= previous day’s road temperature profile. 


** T,= road surface temperature. 


2. Description of input data 

All three models are run on a 24-hour cycle, 
‘projecting forward the measured noon road surface 
temperature. The model inputs are similar for all three 
models and are composed of six elementary meteoro- 
logical variables: air temperature, dew-point, wind 
speed, total cloud amount, cloud type and precipitation. 
Cloud type is assigned | for low, 2 for medium and 3 for 
high cloud. Precipitation is represented with 0 for no 
rain, | for light rain and 2 for moderate or heavy 
continuous rain. For the Met. Office model, low-cloud 
amount is also needed. For comparison purposes the 
mean values over a 3-hourly interval are used. Besides 
the meteorological variables, road surface and depth 
(0.30-0.45 m) temperatures at noon are required by all 
three models. 

To give an objective comparison, the models are run 
against each other in two different ways. Firstly, actual 
data measured at the roadside are used for the input data 
to produce a retrospective prognosis (RSP) — sometimes 
called a perfect prognosis. This enables the internal 
errors of the model schemes (physical and mathematical 
approximations) to be compared, eliminating errors 
related to forecast input. The input for the RSP consists 
of 95 days of data from the winter of 1988/89. Air 
temperature, dew-point and wind speed are taken from 
the measurements of the instruments mounted in and on 
roadside screens at a non-standard height of 2.0 m (for 
temperatures) and 3.0 m (for wind speed). Cloud data 
and rainfall intensity and duration are derived from the 
observations of the nearby weather station at the 
University of Birmingham. 

Secondly, to assess the size of errors due to model 
interior shortcomings plus forecast input errors, the 
models are compared using real-time input data — i.e. 
real-time prediction (RTP). The real-time input was 
produced by Birmingham Weather Centre at Elmdon 


for Hereford and Worcester County Council as part of 
an ‘Open Road’ commercial contract. For the purpose 
of this comparison the days selected are those when the 
forecast input was issued around noon and not updated 
later, which gives a data set of 65 days. 

Each model is run at an initial time of 1200 UTC and 
compared with hourly measurements of road surface 
temperature obtained from a road surface sensor in the 
fast lane of the motorway. The accuracy of this sensor is 
approximately + 0.5 °C, and the calibration was checked 
twice during the winter period considered. 


3. Comparative results for RSP 


3.1 Diurnal and overall analysis of RSP 
The hourly differences of predicted and measured 
road surface temperatures are first calculated for the 
Thornes, Met. Office and Icebreak models with the 
actual observed inputs. The biases and SDs for each 
hour are shown in Figs 1-3. The results can be expressed 
as follows: 
(a) The Met. Office model has a consistent cool bias 
that is more than —0.9 °C for every hour except 1300 
UTC and is as low as —1.5°C in the evening and 
morning. The Thornes predictions are less cool, its 
bias reduces in the late evening and early morning 
(2100-0700 UTC), but increases notably in the late 
morning. The Icebreak model has a relatively small 
bias that is the closest of the three to zero. 
(b) For hourly SDs, all three models have smaller 
values at night. This can be explained by a simpler 
parametrization procedure at night as the main 
external forcing term — solar radiation (which is 
difficult to model due to the unknown shadowing 
effect of broken cloud) — is excluded, and the models 
can cope with the surface energy balance better. The 
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Thornes model has the largest SD, more than | °C 
during the day and night, whilst the Met. Office 
model and the Icebreak model have less than 1 °C 
SDs (as low as 0.7 °C for the Met. Office model and 
0.6 °C for the Icebreak model) at night. 


The overall bias, SD and r.m.s. error of model 
predictions were calculated by averaging the hourly 
errors and are shown in Table III. For the purpose of 
comparison, it should be noted that a bias expresses the 
‘mean status of prediction errors, such that a negative 
bias means that the model predictions are too cold. The 
SD indicates the range of variance of the errors and the 
r.m.s. error is referred to as a measure of the total errors. 
A good model is expected to have a small absolute bias 
and a small SD or small r.m.s. From Table III and 
Figs 1-3 it can be seen that, in the prediction of overall 
temperatures, the Icebreak model is of the highest 
accuracy. The Thornes model has a larger standard 
deviation and thus differs more from actual values. The 
Met. Office model tends to be much cooler than reality 
and has the largest r.m.s. value. For winter maintenance 
applications a small negative bias is preferable — erring 
on the side of caution, i.e. some de-icing chemicals 
wasted but fewer potential accidents (Thornes 1989). 
Obviously, small values of SD and r.m.s. error are also 
required. 


3.2 Minima and time of freezing of RSP 

One of the important aspects of model comparison is 
also to look at the minimum temperature prediction and 
the time of start and duration of freezing. The daily 
predicted minimum road surface temperature for each 
model has been compared with the observed minima to 
yield the mean of the difference of minimum temperature, 
r.m.s., SD and frequency distribution of the difference. 
When both predicted and measured temperatures fall 
to, or below, 0 °C (7,< 0 °C), the time and duration are 
compared and their differences calculated. Table IV 
gives a comparison for the three models in the 
differences of the minima, start time and duration of 
freezing. 

From Table IV it is interesting to note that in the 
prediction of minimum surface temperature, the Thornes 
model has a small bias with a large SD and the Met. 
Office model has a smaller SD with a much larger bias. 
The Icebreak model has both the smallest bias and SD. 
In comparing freezing time and duration, the prediction 
of the Thornes model is 37 minutes later than the actual 
freezing time and 19 minutes ahead of actual surface 
temperature rising above 0 °C. The Met. Office model 
tends to issue an icing warning almost an hour early, and 
the warning lasts 2 hours longer than the actual period 
of freezing. Compared with the Thornes and Met. Office 
models, the Icebreak model has the smallest errors in 
predicting freezing time (30 minutes early). 

The frequency distributions of the difference between 
actual and predicted minimum temperatures for the 
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Figure 1. Hourly error analysis of retrospective prognosis for the 
Icebreak model (Chapman’s Hill, winter 1988/89, 95 days). 
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Figure 2. As Fig. | but for the Thornes model. 
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Figure 3. As Fig. | but for the Met. Office model. 
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Table lll. Overall errors of retrospective predictions for 
model comparison (Chapman’s Hill, winter 1988/89, 95 days) 
calculated from the average hourly errors 








Icebreak Thornes Met. Office 
No. of hours* 2164 2164 2164 
Bias —0.29 —0.61 —1.14 
SD 0.956 1.284 1.048 
r.m.s. 1.000 1.418 1.543 





* Note that 2164 hours of data were available for analysis out 
of 2280 (94.9%) owing to data loss from road sensor. 
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Table lV. Differences of predicted and measured road 
surface temperature minima (°C), freezing start time and 
duration (min.) for model comparison (Chapman’s Hill, 
winter 1988/89, 95 days) 








Icebreak Thornes Met. Office 

Minima: 

Number 95 95 95 

Bias —$ii —0.17 —0.92 

SD 0.757 1.022 0.786 
Freezing: 

Number 19 16 19 

Start time* 30 =a 57 

Duration** 37 —19 123 





* A positive number indicates that the road surface temper- 


ature was predicted to fall to 0 °C before it actually did. 
** A positive number indicates that the duration of time at or 
below 0 °C was predicted to be too long. 


models are shown in Figs 4-6. For an accuracy in 
absolute error in the prediction of minimum road 
surface temperature of + 0.5 °C, 57.9% of the Icebreak 
model, 42.1% of the Thornes model and only 25.3% of 
the Met. Office model are found. Most of the absolute 
prediction errors of the Icebreak model (87.3%) and 
Thornes model (76.3%) are less than 1.0 °C;for the Met. 
Office model predictions 55.5% are of this accuracy. 

An analysis of ‘frost’ (defined simply as road surface 
temperature falling to 0°C or below) or ‘no frost’ 
forecasts against ‘frost’ or ‘no frost’ occurring (forecast / 
actual) has also been made and the results are given in 
Table V. As classified by Thornes (1989), the potential 
consequences of an erroneous forecast are: 

Type 1 error (T1): Forecast ‘no frost’, actual ‘frost’ 

— potential for accidents, 

Type 2 error (T2): Forecast ‘frost’, actual ‘no frost’ 

— potential for wasting de-icing chemicals. 


Table V shows that for the 95 nights (with 21 frosts), 
T1 is the highest for the Thornes model, while T2 is the 
highest for the Met. Office model. 


3.3 Maximum temperature of RSP 

An accurate prediction of maximum road surface 
temperature (usually afternoon temperatures) is another 
important aspect of a successful model. It acts as to 
convince highway engineers of the model accuracy, for 
the engineers compare model prediction with sensor 
measurements in the afternoon, and may believe that an 
incorrect prediction of maxima will lead to a wrong 
prediction of minima in the night-time. 

The analysis of model retrospective prognoses for the 
maximum temperature is given in Table VI. Both the 
Thornes model and the Met. Office model give much 
cooler bias than the Icebreak model, and the Thornes 
model has the largest SD. It is seen from the table that 
the Icebreak model has significantly improved the 
prediction of maximum temperature. 
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Figure 4. Distribution of the difference of forecast minus actual 


minimum temperature in retrospective use of the Icebreak model 
(Chapman’s Hill, winter 1988/89, 95 days). 





<*-3.0 
-3.0/-2.5 
-2.5/-2.0 
-2.0/-1.5 
-1.5/-1.0 
-1.0/-0.5 
-0.5/ 0.0 
0.0/ 0.5 
0.5/ 1.0 
1.0/ 1.5 
1.5/ 2.0 
2.0/ 2.5 
2.5/ 3.0 
> 3.0 


Temperature forecast error (°C) 











° 
= 
°o 


20 30 40 50 
Percentage 


Figure 5. As Fig. 4 but for the Thornes model. 
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Figure6. As Fig. 4 but for the Met. Office model. 


4. Comparative results using ‘realistic’ 
Open Road input data 

To see what the potential accuracy of model 
prediction is in view of the fact that perfect predicted 
data are never available, the three models were run using 
the ‘realistic’ input ‘data used by the Birmingham 
Weather Centre for the nights concerned. This gives an 
estimate of the scale of model prediction errors in 
routine operation. Note that although only 65 days were 
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Table V. Model comparison of frost forecast accuracy with actual roadside 
input (Chapman’s Hill, winter 1988/89, 95 days). (See text for explanation of 








Tl and T2) 
Forecast/ actual Icebreak Thornes Met. Office 
No. % No. % No. % 

No frost/no frost 73/74 98.6 74/74 100.0 69/74 93.2 
Frost/no frost T2 1/74 1.4 0/74 0.0 5/74 68 
Frost/ frost 19/21 90.5 16/21 76.2 19/21 90.5 
No frost/frost T1 2/21 9.5 5/21 ~ 23.8 2/21 9.5 
Overall 92/95 96.8 90/95 94.7 88/95 92.6 








Table VI. Model comparison for road surface temperature 
maxima prediction (Chapman’s Hill, winter 1988/89, 95 days) 








Icebreak Thornes Met. Office 
Number 95 , 95 95 
Bias —0.56 —{.22 —1. 
SD 1.597 1.972 1.571 





selected for this analysis, the results for the retrospective 
model runs given in Table VII are very similar to those 
discussed in section 3. 


4.1 Errors in predicted input 

There is a need to investigate how accurate was the 
forecast input data used in the real-time prediction. The 
mean and maximum differences, SD and r.m.s. of the 
differences between real-time input (forecast) and 
retrospective input (actual) for screen air temperature 
and dew-point, wind speed, cloud amount, cloud type 
and rainfall for the 65 days are given in Table VIII. The 


differences are calculated from the values pertaining at 
3-hourly intervals. It should be noted that there are 
height differences between actual (2.0 m) and forecast 
(1.2m) air temperatures and dew-points, and actual 
(3.0 m) and forecast (10.0 m) wind speeds. 

It is seen from the table that the most reliable 
variables predicted in real-time input seem to be air 
temperature and cloud type. All the others have a large 
difference. The difference in wind speed may be partly 
explained by the difference in the height where wind 
speed was derived, but this can only be surmised. 


4.2 Comparison of model potential 
accuracy 

The purpose of this exercise is to compare the size of 
errors due to model interior shortcomings (as discussed 
in section 3 for perfect predictions) with those due to 
forecast input errors. The comparison of retrospective 
prognoses and real-time predictions (of the same 65 





Table Vil. 


Model comparison of retrospective (RSP) and ‘real-time’ (RTP) 


runs (Chapman’s Hill, winter 1988/89, 95 days) 








Icebreak Thornes Met. Office 
RSP RTP RSP RTP RSP RTP 
No. of hours* 1480 1480 1480 1480 1480 1480 
Overall: 
Bias (°C) —0.33 —0.10 —0.61 —0.17 ~1.15 -§42 
SD (°C) 0.75 1.14 0.97 1.12 0.81 1.02 
r.m.s. (°C) 0.93 1.43 1.31 1.42 1.45 1.28 
Maximum: 
Bias (°C) —0.56 —0.16 ~—1.0 -6.73 —1.13 -0.78 
SD (°C) 1.58 1.63 1.96 1.87 160 1.52 
Minimum: 
Bias (°C) —0.17 —0.38 —0.19 —0.11 —0.94 —0.26 
SD (°C) 0.69 1.64 1.10 1.48 0.85 1.32 
Time of freezing: 
Freezing (nights) 17 17 15 15 17 16 
Start (minutes) 37 64 —28 52 46 —49 
Duration (minutes) 52 56 {2 68 103 0 





* 1480 hours out of 1560 possible (95% availability) 
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Table Vill. Comparison of meteorological variables between 
real-time input and actual roadside input (Chapman’s Hill, 
winter 1988/89, 65 days) 





Realistic with actual 





Mean SD Maximum 

error error 
Air temperature (°C) 0.07 1.534 5.50 
Dew-point (°C) 0.44 1.710 5.90 
Wind speed (knots) 2.75 3.860 17.00 
Cloud amount* (oktas) —0.43 2.236 8.00 
Cloud type* (1—3)** —0.19 0.748 2.00 
Rainfall* (0—2)** —0.53 0.794 1.00 





* Although these are discrete rather than continuous 


variables, the analysis has still been done for comparison. 
** See text for explanation. 


days) for the three models is shown in Table VII. Note 
that the days chosen for this analysis were the 65 
occasions when no update was issued, and hence the 
original noon forecast was considered good, and 
weather conditions did not change significantly from the 
forecast. 

In terms of overall RTP errors, all three models have a 
negative bias and the Met. Office model is the coldest. 
Their SDs and r.m.s. errors are not much different. For 
maximum temperature, the Icebreak model is marginally 
colder than actual maximum (—0.16°C) with SD of 
1.63 °C. The Thornes model and Met. Office model are 
cooler. In the real-time prediction of minimum temper- 
ature, the Icebreak model and Thornes model tend to 
issue a freezing warning nearly an hour earlier than 
actual freezing, while the Met. Office model has a lag of 
49 minutes in predicting the start of freezing. The Met. 
Office model has an error of zero in the prediction of 
duration of freezing conditions. 

The results of forecast/ actual frost analysis for the 
realistic prediction is listed in Table IX. It can be seen 
from the table that the overall accuracy of the three 
model is similar. 

Comparing the ten left columns with the ten right 
columns in Table VII for each model, the Icebreak 
(7/10) and Thornes (6/ 10) models show an improvement 
with retrospective actual inputs, while the Met. Office 
(3/10) model shows little improvement. This suggests 
that the forecasters at Birmingham Weather Centre 
were varying the inputs to produce better results, as 
discussed further below. For the Icebreak model, the 
overall mean bias and maximum temperature bias of 
RSP are a little cooler than those of RTP; significant 
improvements of RSP are seen in the reduction of r.m.s. 
and SD values of overall and minimum temperatures. It 
is clear that with actual input, the model has improved 
its overall and minimum predictions, and freezing time. 
But it shows no improvement in the prediction of 
maximum temperature and freezing time over real-time 
running. 





Compared with real-time predictions, the Thornes 
model shows an improvement in the prediction of 
overall and minimum temperature and freezing time 
with actual roadside input, while the Met. Office model 
has cooler predictions using perfect prognosis. 

Looking at the overall errors, it can be seen that the 
SD of the Icebreak model is reduced from 1.14 to 0.75 
using the perfect prognosis, a reduction of 34%. For the 
Thornes model the percentage reduction is about 13%, 
and the Met. Office model about 21%. The Met. Office 
model shows little potential for improvement using 
accurate or actual input in terms of both bias and r.m.s. 
error however. It is notable that the input from the 
weather forecasters significantly improved the output of 
the Met. Office model. 


5. Discussion and conclusion 

Although it is indicated in the previous section that 
the Icebreak model and Thornes model achieve 
improvements with actual roadside input, perfect 
prognosis for the Met. Office model is not superior to 
realistic prediction — this requires further investigation. 

The first question is whether the difference of 
temperature and wind speed height between actual 
(2.0 and 3.0 m) and forecast (1.2 and 10.0 m) input can 
be considered as a significant contributing factor to 
the abnormal results of the Met. Office model. The 
difference of temperature height is about 0.8 m. It is 
reasonable to believe that this difference of temperature 
height can be ignored because the near-ground layer 
over the road is such a well-mixed layer that no big 
difference of vertical temperature should exist within so 
small a vertical distance, except on exceptional nights. 

The difference in wind speed height is relatively large. 
Because of the inaccuracy of interpolation for a wind 
speed profile, roadside wind speed was not extrapolated 
to 10 m height in retrospective prognoses. By the same 
consideration of a well-mixed layer, however, the 
difference of wind speed due to the difference in height 
is not expected to be significant to the results, as sensi- 
tivity analysis shows that wind speed has a small 
influence on model output (Shao 1990). Therefore it can 
be concluded that the results in the model comparison 
are reliable. 

A probable explanation for the abnormal results of 
the Met. Office model is that the model has been used 
operationally for 2 years and the local forecasters at 
Birmingham Weather Centre have become familiar with 
optimizing model output by varying model input. 

Farmer and Tonkinson (1989) have undertaken a 
verification of the Met. Office model by comparing the 
results of perfect prognoses against actual road surface 
sensor (fast lane Chapman’s Hill) measurements in the 
winter of 1988/89. Their perfect prognoses were made 
using actual roadside observations and airfield observ- 
ation from Birmingham Airport. Their results are 
compatible with those derived here. The results all show 
a cool bias of the Met. Office model predictions. 
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Table IX. Model comparison of frost forecast accuracy with realistic Open 
Road input (Chapman’s Hill, winter 1988/89, 65 days). (See text for explanation 


of Tl and T2) 





Forecast/ actual 








Icebreak Thornes Met. Office 
No. % No. % No. % 
No frost/no frost 45/47 95.7 47/47 100.0 45/47 95.7 
Frost/no frost T2 2/47 4.3 0/47 0.0 2/47 43 
Frost/ frost 17/18 94.4 15/18 83.3 16/18 88.9 
No frost/frost T1 1/18 5.6 3/18 16.7 2/18 11.1 
Overall 62/65 95.4 62/65 95.4 61/65 93.8 
References 


By comparing bias, SD and r.m.s. error in a retrospect- 
ive analysis, we believe that we have shown that the 
Icebreak model provides improved accuracy in predict- 
ions of maximum and minimum road surface temper- 
ature prediction, such as in the prediction of maximum, 
and frost frequency. The predictions of the Thornes 
model have a larger variance, and the predictions of the 
Met. Office model are too cold. 

The predictions of the Icebreak model can be further 
improved with more accurate input, such as for wind 
speed and cloud. The Thornes model is of smaller 
potential to improve its predictions, and the Met. Office 
model has no such potential unless the negative bias is 
removed and the forecasters do not subjectively alter the 
inputs. 

The results presented here are restricted by the use of a 
single test site, and research is now continuing for a 
selection of sites across the United Kingdom. 
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Satellite and radar photographs — 11 January 1991 at 0300 


and 0316 UTC 


On 11 January 1991 a frontal band of cloud (F—F-F in 
Fig. 1) was slow-moving, lying across northern France 
and the south-east coast of England. The eastern 
Atlantic was awash with vortices. In particular, cold air 
comma clouds (labelled A, B, C, D in Fig. 1), revealing 
enhanced convection, were running north-east just 


a 


Figure 1. 


north of the frontal band. These brought thunderstorms 
with hail to southern counties of England and Wales 
during the day. 

Comma A had induced a shallow wave on the front, at 
0300 UTC near Luxembourg (see Fig. | inset). Comma 
B had a well-developed head of cirrus outflow and the 


ve 


Photograph by courtesy of University of Dundee 


NOAA-1I1 Channel 4 (infra-red) satellite picture at 0316 UTC on 11 January 1991. F—F-F is a slow-moving frontal band of cloud. 


A, B, C, D are four of the many vortices covering the Atlantic. The surface analysis for 0300 UTC is inset. 
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convection in the tail of the comma can be clearly seen 
beneath the frontal cloud to the west of Brest in Fig. 2. 

The COST-73 satellite/ radar composite for 0300 UTC 
(Fig. 3) shows the line of convective precipitation in the 
tail of the comma. Unfortunately the picture is 
complicated by previous storms to the north of Brittany. 
At 0600 UTC (Fig. 4) the radar picture is clearer, with 
the line of showers and thunderstorms extending well 
south over Brittany beneath the frontal cloud. 

Radiosonde ascents show the air over England and 
Wales was unstable to around 400 mb — Larkhill 
we, 5°W 


953 (Stn No.) 
vy 00 (UTC 


Convection 





i 


(03743) is typical of this air (Fig. 5). Convective tops 
over Brittany were limited to 600 mb. 

The orientation of certain cloud patterns is determined 
by the windflow relative to a moving system, rather than 
the actual winds reported, so it is useful to examine the 
relative flow to clarify the dynamical interpretation of 
the imagery. The relative wind is found by vectorially 
subtracting the system velocity (for comma B this was 
245° 35 kn) from the actual winds. Plotted winds in 
Fig. 2 show the relative winds at 400 mb from radiosonde 
stations, plotted at appropriate positions relative to the 








Photograph by courtesy of University of Dundee 


Figure 2. As Fig. | but Channel 3. For explanation of wind plots and insert see text. 
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Figure 3. COST-73 satellite and radar image for 0300 UTC on 
11 January 1991. Light blue represents cloud tops with temperatures 
between —15 and —45 °C and dark blue < —45 °C. Rainfall intensities 
(mm h'') are shown as follows: pink < 1, green 1-3, yellow 3-10, red 
10-30. Coastlines, national boundaries and the limits of radar 
_ Coverage are shown in black. 





Figure 4. As Fig. 3 but for 0600 UTC. 


comma head at 0300 UTC. In general these show the 
cirrus outflow being carried to the north-east. However, 
Aberporth (which was near the centre of the comma 
head at 0700 UTC) has a relative wind of 135° 24 kn. 
This agrees well with the anticyclonically curved 
striations in the cirrus outflow. Although layered cloud 
is also present due to mass ascent within it, the comma 
head as seen by the satellite is mostly cirrus outflow from 
cumulonimbus clouds feeding into its southern flank. 
This is confirmed by the precipitation shown in the 
radar pictures. 

Another notable feature was that individual storms in 
the Channel could be tracked for over 5 hours. To 
achieve this life-span, new cells were building into each 








225/98 

















220/52 
225/61 
225/53 
di 
i 
di 
230/40 
\ 7245/44) 
, 7 _540/50—4 
235/41 
\ \ \ é 














~% 225/36 
if 4 
a \ A¢ | 268/27 
Z 180/14 
Z\ 
—10° 0°c 10° 20° 


Figure 5. Tephigram for Larkhill (03743) at 0700 UTC on 
11 January 1991. 


























Figure 6. Hodograph for Larkhill (03743) at 0700 UTC on 
11 January 1991. Winds shown are relative to the storm motion. 


storm, and this requires directional shear in the winds 
(see Bennetts et al.*). Fig. 6 shows the Larkhill winds 
relative to a storm moving from 245° at 40 kn. With sea 
temperatures of 10 or 11 °C, convection would be from 
the surface. The low-level inflow is from the south-east 
quadrant, the cells lean to the north, and outflow is to 
the north-east. This can be inferred in some of the 
individual storms near the Larkhill plot in Fig. 2. 
Using isentropic analysis and relative winds, a 
tentative conceptual model of the comma cloud is given 
in Fig. 2 (inset) showing the ascending conveyor belt. 
J.R. Grant 





* Bennetts, D.A., McCallum, E. and Grant, J.R.; Cumulonimbus 


clouds: an introductory review. Meteorol Mag, 115, 1986, 
242-256. 
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